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SUMMARY 


Approximate solutions of the differential equations of continuum 
isentropic gas dynamics may be obtained for flows in the equatorial 
plane of a binary star atmosphere if the analogy between gas dynamics 
and shallow water flows is employed to define pertinent experiments in a 
rotating water tank. The magnitudes of the local inclination of the 
tank's bottom and of the resulting local gravitational force at the 
corresponding point in the binary star atmosphere are proportional. 

The analogy requires a constant polytropic coefficient of two. This 
application of the analogy rests on its extension presented here to 
rotating systems with external force fields. A theoretical evaluation 
of the analogy, for the case of the earth's atmosphere under hydrostatic 
conditions with polytropic coefficient equal to two, yields the correct 
result for the maximum height of the atmosphere where density reaches 
zero. This study is motivated by an application to the early stages 
after a hypothetical fission of the earth and the moon. 


tThe first author is consultant to DVL and to the Goddard Space Flight 
Center of the National Aeronautics and Space Administration. 
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LIST OF SYMBOLS 


A 

B 

c 

Pg 

F. 

jw 

r 


w 


I! 


max 


M 

P 

P, 

r 

r 


oo 


T _ __ 

u,v,w 

x,y,z 

r 

e 

\ 


H 

X 

w 

li 

p 

I 

1 

uo 


characteristic constant defined in (26) 
characteristic constant defined in (32) 
specific heat at constant pressure 

surface function defined in connection with (10) 

FROUDE -number 

gravitational acceleration 

depth of water rotating as a solid body, defined in (13) 

maximum altitude of atmosphere 

mass of celestial body 
pressure 

atmospheric pressure at the free surface of water 

distance from center of rotation 
reference length 

gas constant 

time 

absolute temperature 

velocity components parallel to the coordinate axes 

Cartesian coordinates 

universal gravitational constant 

small dimensionless length, defined in connection with (34) 
disturbance amplitude of water depth 

polytropic exponent of a gas _ 
wave length corresponding to T[^ 
mass ratio 
density 

nondimensional gravitational potential 
gravitational potential 
constant angular velocity 


Subscripts 

reference quantities or quan- 
tities taken at the surface 
of the earth 

bottom of container of water 
free surface of water 
free surface of water in case 
of rotation of water as a 
solid body 


All dimensional quantities are marked by a bar above the letters. 


i gas and water o 

g gas 

w water 

a, 3 first and second celes- 1 

tial body, respectively 2 

II 


1 


INTRODUCTION 


According to observation, a binary star system may possess a 
common atmosphere which envelops both stars. The motion of a particle 
in such an atmosphere may be treated by use of the restricted three -body- 
problem in celestial mechanics. This approach has been employed in 
several publications, e.g., Kopal (1956, 1957) and Mrs. Gould (1957* 

1959)- (See part E of this paper.) According to Prendergast (i960), 
the mean free path in these atmospheres is of the order of magnitude of 1 
- 10 km, i.e., this characteristic length is so much smaller than either 
the separation of the stars or the radius of one of them that collisions 
between particles in the atmosphere have to be taken into account. 
Therefore, Prendergast (i960) studied stationary solutions of the dif- 
ferential equations for continuum isentropic gas flows in the equatorial 
plane of the binary system by neglecting (a) the pressure gradient and 
(b) the velocity component normal to the "Lagrangian surfaces". Sub- 
sequently, Huang (1965) employed a stead-state celestial mechanics 
approach which neglects the pressure gradient but takes into account 
collisions between particles by invoking the statistical properties of 
the Jacobian constants of the colliding particles. 

The theory presented in this paper was developed because of a 
desire to study gaseous motions in the vicinity of earth and moon during 
the early stages after a fission of their hypothetical parent body. By 
such motions, mass, energy, and angular momentum would be redistributed 
and possibly exported from this binary system. 

The relatively large mean free path quoted above may be assumed 
to represent an average value for the entire atmosphere enveloping the 
stars revolving around each other. A rigorous hydrodynamic treatment 
of the atmosphere including the pressure gradient, therefore, seems 
desirable at least for the regions of maximum density within the atmos- 
phere. In addition, the observed eruption phenomena in binary systems, 
e.g., Huang (1965), render desirable a non-stationary treatment of the 
atmosphere. If these atmospheric motions always are 

symmetrical with respect to the equatorial plane, the differential equa- 
tions of motion and energy for three-dimensional isentropic gas flows 
reduce regorously to their two-dimensional versions in this plane; in 
the continuity equation, however, one term has to be neglected in this 
two-dimensional transient approach. Since Prendergast (i960) points out 
that the total mass of the atmosphere is concentrated about the equatorial 
plane, a study of flows in this plane appears to be relatively unimportant 
a restriction of generality which, in addition, governs the preceding 
pertinent publications mentioned above. Even in the equatorial plane, 
the time -dependent nonlinear system of differential equations with three 
independent variables precludes any analytic solutions. Numerical solu- 
tions reveal general characteristics only after a very large number of 
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individual cases have been computed and evaluted. Because of this 
situation, the analogy between gas dynamics and shallow water motions 
is employed in the following to furnish the theoretical background for 
a water-tank-experiment simulating transient flows in the equatorial 
plane of the binary star atmosphere.* This requires a discussion of 
shallow water flews (part B) and an extension of the validity of this 
analogy to flows in rotating systems which are subjected to external 
force fields (part C). A presentation of an analogy such as in this 
paper should include some typical and sufficiently general results. 
Unfortunately, means are not available to the authors to carry out 
these experiments for the time being. Since this situation will not 
change in the foreseeable future, it was decided to publish the paper 
now. As a restricted test of validity and accuracy of the method 
presented here, the hydrostatic density distribution in the earth r s 
atmosphere is obtained as a special application in part D. 

SHALLOW WATER FLOWS - PART B 


A container, schematically shown in Figure 1 is filled to a 
certain height with water. At time t =0, this water tank is rotating 
about the z ^.axis with constant angular velocity The acceleration 

of gravity W g is assumed to act in the negative z^-direction. A cylin- 
der with radius r is mounted at the bottom of the container; its axis 
coincides with the W zL-axis. Along the wetted surfaces of the cylinder 
and/or the bottom, sources or sinks may add water to the system or 
drain it off according to some law, which may vary locally and/or in 
time. Thus the total amount of water in the contained need not neces- 
sarily be constant. The shape of the bottom will be treated later on. 

The subscripts w and g denote the water and the gas, respectively; 
the subscript i stands for both w and g. Dimensional properties are 
marked by an upper bar. The angular velocities nJ^, the lengths r . , and 
the densities p" . will be used as constant reference quantities, rlon- 
dimensional Cartesian coordinates x^,y^,z^, velocities u^,v^,w_^, times t^, 
pressures p. and densities p. are introduced by relations such as the 
following : 


x. 

1 



u. 

1 



t. = 5>. t. 

i li 


p. = ^ _ 2— 2 

i p .0). r . 
*01 1 oi 


p ± 



with p = p* 

r w ow 


const, and TT = p (r ) 
og r g og 


*For general information about this analogy, e.g., Courant/Friedrichs 
(19^8) and Wehausen/Laitone (i960) are recommended. 


( 1 ) 
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'If. = r . <£>. 2 /g\ = const.. 
J 1 Ol 1 ' & x 


( 2 ) 


a FROUDE -number ^ r *. is introduced as the ratio of the centrifugal accelera- 
tion ah 2 at the distance r . from the axis fo rotation to the gravita- 
tional acceleration g. . The symbol D/Dt denotes the substantial deriva- 
tive with respect to £ime t . According w to Lamb (1932, p.4 and p*3l8), the 
equations of motion of the water flow read as follows in a system rotating 
at constant angular velocity ui : 


Du 

=rr ~ “ 2v - x 
Dt w w 

w 



dx 


w 


Dv 
_ w 

Dt 

w 


+ 2u 
w 


- y 


w 


dp 


w 


ay. 


w 


(3) 


and 


Dw 


w 


3P, 


w 


Dt 


w 


dz_ 


w 




(4) 


Here (4) was obtained from 


Dw 


w 


w Dt 


w 


dp 

w - 

dlf ®w 

w 


(5) 


using relations (l). The left hand sides of (3) - (5) express the mass- 
fixed accelerations in a system at rest provided the quantities of the 
rotating system are used. According to Goldstein (i960, p. 12), the 
equation of continuity is independent of the chosen coordinate system, 
since it simply states the constancy of a scalar quantity. It thus reads 
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! 


( 6 ) 


Su 

w 



+ 


5 v 


w 


s y. 


w 


dw 


w 


dz 


w 


= 0 . 


For a) - 0 , the basic simplifying assumptions of the shallow 
water theory are, see e.g., Wehausen and Laitone (i960, p. 667): 

(I) irrespective of flows in the water, the pressure distribu- 
tion in a vertical column of water is the same as in hydro- 
statics and 

(II) the horizontal velocity components u and v are independent 
of z , (considering the no -slip condition a^ the bottom, 
this assumption can be satisfied only approximately and 

for short times . ) 

In addition to assuming the validity of (i) and (.II) for u) ^ 0 , 
here also the following assumption has to be made: 

(III) the hydrostatic pressure distribution holds true even in 
the presence of motions relative to the rotating frame of 
reference . 

If the fluid rotates as a solid body, its free surface (subscript II) 
is given by 


zL-Cr ) = const. + r 2 /2g with r 2 = x 2 + y 2 
II V w' w ^w ' & w w w J v 


( 7 ) 


Should any flows relative to the rotating frame of reference ( "distur- 
bances") be superimposed to this solid body rotation, the free surface is 
given by 


z w Z 2v/ Xw *^w , ^w^ 'V ^w'^w'^w^ + Z II <V*w> 
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If the vertical accelerations Dw /Dt caused by these disturbances are 
negligible due to a sufficiently sraaYl water depth, then ( 5 ) is satis- 
fied by the hydrostatic distribution specified in assumption (III) above. 


P w ?w ®w [ z II^ x w’ y w^ + ^^w'^w’SP z 


w 


+ p. 


oo J 


where p is the atmospheric pressure acting upon the free surface. 
Eq. (8)°8as been derived in Lamb (1932, p. 3l8) for rotating shallow 
water flows. A nondimensional form of (8) is given by 


p w = 


'OO 



z ) + (z 
1W 1W 



The functions F - z - z (x ,y ) = 0 and F - z - z (x ,y ,t ) 

1W V 1 W V 2W w 2w w w w 

-- 0 represent the wetted bottom and the free water surface, respectively, 
see Figure 1. Therefore, 


DF. /Dt = 0 
jw' w 


j = 1, 2 , 


according to Lamb (1932, p. 7)> since these surfaces are always composed 
of the same fluid particles and thus 



Because of assumption (II) and eq. (ll), integration of the continuity- 
equation (6) with respect to in the limits Zi w and gives 


5 ( Zptj - z iw ) 

dt 

w 




(12) 


For tu =0, equations (3) and (12) represent the shallow water 
theory which is a first order approximation by expanding in powers of a 
parameter £ , defined as the ratio of the squares of the characteristic 
horizontal and vertical scales of length, Wehausen and Laitone (1960,__ 
p. 466 and 667). Subsequently, the wave length X and the amplitude of 
the disturbance are employed for these scales . According to Laitone 
(1962), a second order approximation yields in case of X w -* 00 and uS = 0 


9 P. 


w 


dz 


w 


^w®w 




(15) 


where H = z TT 
w II 


- z 


lw 


Correspondingly, Laitone (1962) shows for 0 and <« w - 0 


dp 

- p g | 1 

dz w w [ 

w 



n2 

\ 

J 



and 


(14) 


du 

V 

az 06 

V 


^W H W 


TT H 


W 


W 


y 


r 

w 
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Relations (13) and (14) show for 35 =0 that the following quantities 

have to be sufficiently small in order to allow an inviscid water flow 
to be approximated by the shallow water theory: 

(a) the hydrostatic depth H^(x^y^) = z^ - of the water, 

(b) the ratio T^/H^ of the disturbance amplitude T|^. to H^, and 

( c ) the ratio of H /X 

v ' w w 

It is assumed here that conditions (a), (b), and (c) also govern the 
validity of the shallow water theory in the rotating water tank provided 

both 35 and I grad z (x ,y ) I are sufficiently small, 
w i w w w 

THE ANALOGY BETWEEN SHALLOW WATER AND GAS FLOWS - PART C 

By assuming z = w = 0, a two-dimensional transient gas flow 
(subscript g) will^be rilated to the three-dimensional transient shallow 
water flow (subscript w) treated in part B. By use of H (l) a s defined 
in (13) and taken at station x 2 + y 2 - 1, the following functional 
relationships now are introduced between the water flow and gas flows: 


x 

g 


w 


y g = v 


u 

g 


u ; 
w 



(15) 


P 


g 



- z (.x ,y 
iw w w 


>]/ 


H „ (1 > 


(16) 


and 


p --= j" (x ,y ,t ) - z (x ,y )1 /2^ H (l) = 

g L 2w w v w iw J ' w w' ' 


H „U) 




w 


( 17 ) 
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The scales of length and time of the two flows under consideration are 
connected by the quotients r /r and t ft = (5) /to = const., respectively. 
Since the isentropic flow of a chemically homogeneous gas is considered 
here, there exists a functional relation p = p (p ) between pressure and 
density. The special polytropic relation § <x g p Sh with h = 2, which 
follows immediately from ( 16 ) and ( 17 ), can S only Ipproximately hold true, 
since the kinetic theory of gases shows that 1 < h < 1 . 67 - 

Substitution of (15) - (17) into (5) and (12) yields by use of (9) 


Du 

P “ 2v - x = - 

g LDt g g J 


5a. p 8(z iw /3r w> 

dx H g §x ’ 


( 18 ) 


-Dv 

) ■=r^- + 2u - y l = 

g LDt g J g J 


5a_ p 

\ 8 3y g 


(19) 


and 


dp 

dt 


+ 


g 


d 

dx 

g 



V + 





0 . 


( 20 ) 


The gas flow defined by (15) - (17) is two-dimensional, since the third 
coordinate z of the corresponding three-dimensional water flow is needed 

to describe ?he density p of the gas. Because of 

§ 


it does not follow from (6) and (20) that p is constant even though this 
seems to be suggested by relations (l) and s (15). By the argument 
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following equation (5)* the left -hand-sides of (l8) and ( 19 ) are the 
mass-fixed accelerations in the rotating system of reference, and (20) 
is the two-dimensional equation of continuity of a compressible medium 
in this system. The last terms in the right hand sides of both (l8) 
and ( 19 ) may be interpreted as due to a conservative body force, which 
can be derived from a potential z (x , y )j^T • Thus (l8) - (20) 
are the equations of motion and cifftiSuitf (in the continuum regime) 
of two-dimensional gas flows subjected to conservative body forces in 
a system rotating at a constant angular velocity. 

The first law of thermodynamics yields for any mass element of the 
supposedly isentropic gas flow the relation c dT = d? /p between __ 
the specific heat ‘c at constant pressure, tS§ abiolute® temperature T , 
the pressure p , anfi g the density p of the gas. Therefore, the energy® 
following ( 5 ) ® ® 


P 


g 


DT 



% 

DT 

g 


( 21 ) 


Eq. (21) is satisfied by any two-dimensional transient flow provided a 
thermally perfect gas and the polytropic relation p cc p" K with k = const, 
are specified. Thus the relation between water and^gas flows as given 
by ( 15 ) _ (IT) is a complete analogy within the limitations of the assump- 
tions mentioned at various placed above. 

The generalization of the analogy developed above consists of ad- 
mitting (a) conservative body forced which determine the bottom shape 
z w = z w (x , y ) of the water tank and (b) rotating systems. Variable 
bottom configurations have been used in many applications of the analogy 
to problems without body forces in order to admit arbitrary constant 
values of h. During the years 1952 to 1962, this issue initiated a dis- 
cussion in which body forces were not mentioned at all. In agreement 
with the derivation presented above, Laitone (1952, 1953* 1961, 1962) 
and Wehausen and Laitone (i960) pointed out the necessity of z lW (x w , y^) 

= const, in case of problems without body forces which were considered 
solely by them, Loh (1959* 1962) and Bryant (i960, 1962) remarked that 
z iw( X w’ ^ ^ ^ cons "k’ i s admissible in view of the other approximations 
involved in the analogy, such as constant values of k* 

The conditions for the validity of the shallow-water theory mentioned 
at the end of section B also apply to the analogy. Experiments by Laitone 
(1952) confirmed that the analogy yields satisfactory agreement between 
gas dynamics and water flows as long as for tap water with usual surface 
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tension^the water depth is smaller than one quarter inch and the wave 
length \ of the disturbances exceeds 3 inches, see also Laitone (1953, 
1961, 19o2). Under consideration of the surface tension of water, Gupta's 
analysis (1965) yields about the same limits for the validity of the 
analogy carried out in water tanks. 

For sufficiently small values of u^, the wave systems caused by the 
tank's rotation (e.g.. Lamb (1932)) presumably are negligible as compared 
to the wave systems which have been discussed above for uj =0. Since, 
in addition, a transient boundary layer growth at the tanS bottom is due to 
the waves disturbing the solid-body rotation of the water, the relation 
t Ho - t ou following from (1) and (15) prohibits atmospheric motions 
t§ b§ stuSieS in the water tank beyond sufficiently short time intervals 
t after the beginning of the said waves. The required small values of 
dr demand a sufficiently small maximum inclination of the tank bottom. 

This may necessitate the addition of a detergent to the water because of 
the surface tension. 

APPLICATION TO AN ASTROPHYSICAL AND TO A GEOPHYSICAL PROBLEM - PART D 


The right-hand sides of (l8) and (19) were interpreted as representing 
body forces, which result from a potential z lW (x , y ) / f . The gravi- 
tational potential of spherically symmetric celeftia? bodies now will be 
employed to demonstrate the analogy. Figure 2 shows a "binary" system 
of two celestial bodies rotating about their common center of gravity S 
a,t the constant angular velocity $ . The masses of the bodies are M and 
M_, r and r denote the distance §f any point Q(x,y,’z) from the centers 
or th^ two masses, respectively. The gravitational potential in the 
space outside M and Mp then is given by 


Y = 


f M 
a 

r 


a 



( 22 ) 


where F is the universal gravitational constant. In the following, the 
contribution of the atmospheric matter to the gravitational potentials 
will be neglected. The mass ratio h = M./M is introduced and the radius 
of any one of the two celestial bodies is uS'ed as a constant reference 
length r , i.e., r /T and r = r./r . Since T has the units (cm 2 sec' 2 ) 
the expreision r 2<X u5°§ = const, rs employed to obtain nondimensional 
potential ® 
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(23) 


Y = 




where the Frou.de -number 


7 


=r 3 — 2 
r (JO 

og ~ 


T TT 


a 


is the ratio of centrifugal to gravitational forces at the surface of 
the celestial body M . In applications of the analogy, t - r- The 
components of the gradient of this potential now will be SquateS. to the 
last terms in equations (l8) and (19) to obtain an expression for the 
bottom shape in the water tank 

z 

- grad = grad Y • 


With z 

o 


as a constant of integration, it is found that 


z 

0 


-z 

iW 



(24) 

(25) 


g; 


The analogy as expressed by (15) - (17) enables one to study two- 
dimensional atmospheric motions subjected to a gravitational field T(x^, 
by experimental simulation in a water tank with the bottom z -- z (x , 
These atmospheric flows are restricted to the equatorial plane z =0 
where the complete differential equations of motion and energy tike on 
the special two-dimensional forms (l8) - (19) and (21), respectively, 
provided the motion is isentropic and symmetrical with respect to the 
plane z =0. Due to this symmetry, the derivatives with respect to z 


w 


of 


u 




and T vanish identically as functions of x ,y and t 
ii tie flaie z = ®0. Since dw /dz =j= 0 for z -- 0, the^thfee-diminsional 
equation of continuity retains fhis®term dw /d§ which is absent in (20). 
If g s 


(aw /azj 


= 0 
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is known as a function of x ,y and t by use of some information outside 
of the analogy under discusfioS, this S term can be represented by a suit- 
able distribution of sources and sinks in the bottom z = z (x my ) of 
the water tank. The equations of continuity (6) and (?0) tften possess 
corresponding additional source terms. 

Since an experimental verification and test of the accuracy of the 

analogy can not be carried out by the authors for the time being, a 

theoretical check is desirable. Such a test of accuracy is provided by the 

application of the analogy to the earth's atmosphere under conditions of 

a known density distribution, i.e., in case of the hydrostatic atmosphere. 

In this case, the above equations simplify somewhat since y, = 0 and 

F M /r 2 = g , where g is the gravitational acceleration at the surface 
a og & o* & o 

of the earth. Then the quotient / reduces to 

w a 


w 

a 


r u) 1 
ow w 

r ou 
og g 


£• = A , 


(26) 


i.e., the ratio of the centrifugal forces in the container to those at 
the earth's equator. By realizing that r^ -= r^ = r^ because of (15) and 

by taking z = 0 at r =1, equation (25) yields 

lW w 


z 


lW 


U " r ) 
w 


( 27 ) 


The atmosphere enveloping the earth is assumed to rotate as a solid body, 
such that in the analogy z = z IIw and = 0. From (7) it is thus found 


z 


IIw 


= const. 




(28) 


15 


Because of (l6) and (17), 


(29) ^ 

P cr 



J lw 


). 


The equation of state for a thermally perfect gas, p = p R T q , requires 
at the surface of the earth (subscript o) ' " L ° ' 


P R 

(30) 126 = -& 

p °g g Q 

Because of 7 = 

Of 


(31) 

hTv'r =1 = 


r w X 

■with 


(32) 

R T 

B = ~£_^£ 


g r 
0 og 

From 

(27), (28), 

(33) 

Z IIw Z lw 


T 

_2S i_ . 

r ^ 
og "'g 

^ and Z-. = 0 at r 

g Tw w 

2 A B 


and (3l) finally 

= A f2B + |& (r 2 -l) 

L W 


= 1, the expressions (29) and ( 30 ) yield 



is obtained to express the nondimens ional density distribution in the earth ! s 
atmosphere under hydrostatic conditions and for k = 2 . From left to right, 
the terms in the bracket of (33) represent the contributions of the earth's 
surface temperature, of the centrifugal force due to the assumed rigid-body- 
rotation, and of the gravitational force. 


Because of the small ratio of the maximum extension of the atmosphere 


to the earth's radius, it is reasonable to write r 

* r 


- r 


w 


= 1 + e* where second 


and higher powers of e are to be neglected. By realizing that l/(l+e )=l-e+e 
the following relation is obtained to replace (33): 


(34) z IIw - z lW = A [2B - (1 -Xg) «]. 
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Because of (17 ) , the pressure variation in the earth's atmosphere 
according to the analogy is 



Z TT - z 

IIw jW 

(z TT j 

IIw iW 0 


i -r * 

(i 2B ^ ^ • 


( 35 ) 


where the effect of the centrifugal force is negligible even at the 
earth's equator since Jr~ « 3-32 x 10 3 . Because of B = 1.27 x 10 3 

S Q . 

for a surface temperature T = 293 ( K) , the analogy predicts that 
the atmospheric pressure "p ^vanishes at a maximum altitude as given 
by S 


H = r . c = r 
max og max og 


2B 1 - 

1 -T ** 39^ r og 

o 


1 6 km 


It will now be shown that this result for H agrees with the 

max ° 

one of the well-known hydrostatic polytropic atmosphere with k =2, 
where is defined by (a) neglecting rotation, i.e., 3^ =0 and (b) 
by assuming a constant acceleration of gravity. Undfr these condi- 
tions Prandtl/Tietjens (p. 33) show that 


K 

Ja _ f 1 _ H-l h , j H-l 
Po g L H h 0 J 


(36) 


where in the notation of this paper 


h 



r 

°g 


and h = B r 

o og 


(37) 
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Because of e = h/r , Eq. (26) yields in case of h = 2 the relationship 

og 


- y (38) 


which is identical to Eq. (35) for = 0. It should he recalled 

here that the polytropic hydrostatic ^atmospheres defined by condi- 
tions (a) and (b) mentioned above possess a finite maximum altitude 
S' , with the exception of h = 1. 

TT1SLX 

OBSERVATIONS AND GENERAL CONSIDERATIONS WITH RESPECT TO THE ATMOSPHERE 
OF BINARY STARS - PART E 


Prominance activity in the sun is a well known phenomenon. 
Observation of binary stars suggests that the prominence activity 
of the secondary component is stronger on the hemisphere facing the 
primary star. Appropriate initial velocity vectors of these promin- 
ences could result in a gaseous stream toward the primary or escape 
from the binary system, Sahade (i960). 

Within the context of the restricted three -body -problem of celestial 
mechanics, several authors have studied trajectories of ejected 
particles in the equatorial plane. For a particle moving in this 
plance, the Jacobi integral shows the existence of curves of zero 
velocity in a frame rotating with the binary system. These closed 
curves are the "Roche equipotentials" on which the geopotential is 
constant, i.e., the sum of the gravitational potentials of the com- 
ponent stars and the potential of the centrifugal force due to the 
rotation of the system. These Roche equipotentials define admissible 
cross sections of each component star in the equatorial plane . The 
equipotential enclosing the largest area in this plane and only one 
component star is called the Roche limit. For a particle having a 
given value of the Jacobi integral, the corresponding Roche equipoten- 
tial represents a barrier which the particle, of moving inside its 
folds, can never penetrate. The trajectories of particles ejected 
from a component star completely filling its Roche limit take place 
in an area of the equatorial plane whose Roche equipotential boundary 
envelopes both component stars. For this case, trajectories of 
particles have been computed by numerical integrations of the equa- 
tions of motion of the restricted three -body -problem, Kopal (1959)* 

The results could be summarized as follows: (l) particles fall upon 

the companion star, (2) particles fall back upon the parent star, 
and (3) particles escape from the system if they possess sufficiently 
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high initial velocities, generally of the order of several hundred 
kilometers per second. 

According to Sahade (i 960 ), velocities of 'JOO km/sec have been 
observed in the binary system V 444 Cygni. Because of the smaller 
masses, considerably smaller velocities would be needed for the escape 
from the earth-moon system after fission of its hypothetical parent 
body. 


These results of celestial mechanics referred to above do not 
account for interactions between particles. Struve (1957)> therefore, 
has pointed out that it would be important to consider the problem 
at hand also from a hydro dynamical point of view and to allow for 
the effect of radiation pressure and magnetic forces. The theory 
presented in Parts B, D, and D of this paper furnishes a background 
for a water-tank experiment which attempts to simulate gaseous 
streams in the equatorial plane of a binary star system by accounting 
for gravitational, inertial, and hydrodynamic pressure forces. These 
results would pertain to the continuum flow regime which is separated 
from motions due to applications of celestial mechanics by the rarefied 
flow region. The transient water-tank experiment, therefore, would 
furnish a necessary supplement to published results from (a) applica- 
tions of celestial mechanics, Kopal (1959) a nd (b) stationary continuum 
flow theory under rather restrictive assumptions (see part A), 
Prendergast (i 960 ). 

COMMENTS WITH RESPECT TO THE EXPERIMENTAL SIMULATION OF THE FLOW IN 
THE EQUATORIAL PLANE OF THE BINARY STAR SYSTEM - PART F 

In the water tank, the primary and the secondary stars are 
represented by cylinders mounted parallel to the table's z-axis of 
rotation, see Fig. 3. Since the primary star by definition possesses 
the larger mass, the bottom of the tank reaches its smallest z-coor- 
dinate at its intersection with the "primary" cylinder. If the sur- 
faces of the cylinders are sufficiently smooth, the motion of the 
water is very little affected by rotations of the cylinders about 
their own aces of symmetry. Therefore, the wetted portions of these 
surfaces may have to be roughened artificially if this eigen-rotation 
of the cylinders is to influence appreciably the water flow. Even 
though the gravitational fields extend to infinity in the equatorial 
plane, it is sufficient to place a board parallel to the axis of 
rotation of the system on the bottom of this tank, provided the board 
is located beyond the "synchronous satellite orbit" of the system. 

A water sink may be placed at the intersection of the board and the 
bottom. 

Under hydrostatic conditions, the water depth at the circum- 
ference of one of these cylinders may be prescribed arbitrarily. 
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Because of (32) and (33), this depth is proportional to the tempera- 
ture of the atmosphere at the surface of the corresponding component 
star. This temperature should he sufficiently small to ensure the 
gaseous state of the atmosphere since ( 1 6 ), ( 17 ), and the equation of 
state of a thermally perfect gas show that the temperature of the 
atmosphere is proportional to its density in the analogy being em- 
ployed here. The tank's bottom between the cylinders does not have 
to be covered by water continuously in time. As far as mass transfer 
from the secondary to the primary star is concerned, it is sufficient 
to have intermittently water jets emanate from the "secondary" cylinder. 
The mass flow from the secondary to the primary cylinder can be visualized 
by adding dye to the water pouring out of the secondary cylinder. 

CONCLUDING REMARKS - PART G 

The objective of conducting a water-tank -experiment should be a 
qualitative simulation of transient astrophysical effects such as 
the ones mentioned in Part E of this paper. The following details 
militate against quantitative validity of the results: 

(1) the assumption of a constant polytropic coefficient (of 2) 
which is essential to the analogy in every possible version: 

(2) the assumption of isentropic gas flow; 

(3) the treatment of flows in the equatorial plane without regard 
for the interaction with the general three-dimensional flow field 
outside of this plane; this shortcoming manifests itself in 
neglecting the term Sw/dz in the three-dimensional continuity 
equation; 

(4) the omission of radiation pressure, turbulence, viscous friction, 
heat conduction, and other real gas effects; 

(5) experimental difficulties, e.g, such as the ones discussed 
at the ends of part B and part C. 

Even though this seems to be a formidable accumulation of shortcomings 
of the method presented in this paper, the approach may be expected 
to yield more reliable and realistic results in sufficiently dense 
atmospheres than the theoretical methods employed by other authors, 
see Section A. 
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Figure 1: Water Tank. 




Figure 3: Schematic configuration of water tank 

suggested to study the atmosphere around 
the celestial bodies shown in figure 2. 



